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We theoretically model a nuclear-state preparation scheme that increases the coherence time of 
a two-spin qubit in a double quantum dot. The two-electron system is tuned repeatedly across a 
singlet-triplet level-anticrossing with alternating slow and rapid sweeps of an external bias voltage. 
Using a Landau-Zener-Stiickelberg model, we find that in addition to a small nuclear polarization 
that weakly affects the electron spin coherence, the slow sweeps are only partially adiabatic and 
lead to a weak nuclear spin measurement and a nuclear-state narrowing which prolongs the electron 
spin coherence. This resolves some open problems brought up by a recent experiment [f4|. Based 
on our description of the weak measurement, we simulate a system with up to n = 200 nuclear spins 
per dot. Scaling in n indicates a stronger effect for larger n. 
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Introduction — Since the electron spin in quantum dots 
(QDs) has been proposed as a qubit [IJ, much progress 
has been made to develop reliable semiconductor devices, 
mostly with GaAs, in which only a single electron can 
be confined and its spin can be controlled Despite 
the achievement of reliable control, spin decoherence due 
to the hyperfine interaction with the surrounding nuclear 
spins remains a major problem In double QDs, where 
the singlet (S) and one of the triplet (Tq) states can be 
chosen as an effective two-level system to implement a 
qubit 0, H, @]j it has been shown that the probability, 
averaged over many runs of the experiment, to find the 
system in the singlet state at time t > 0, having prepared 
a singlet at time i = 0, is a decaying oscillating function 
(X cos(Jt)exp(-(i/r2*)^) 0,11, with singlet-triplet ex- 
change splitting J and decoherence time 

An important goal in the quest to overcome decoher- 
ence of spin qubits in solid-state devices is to find mecha- 
nisms that allow for an increase in T2. Spin echo has been 
used to reveal T2 ^ fis^ Tj* [8] which sets the scale that 
can be achieved by nuclear-state preparation. In princi- 
ple, it is possible to reduce the nuclear fluctuations, thus 
prolonging , by projecting the nuclear state into (ap- 
proximate) Overhauser eigenstates with either electrical 
or optical means [9|, tlOi] or by polarizing the nuclear spins 
[Tllll^. However, a sizable enhancement of electron spin 
coherence would only be realized for a polarization of 
more than 99% [l^, which so far has not been achieved. 

Reilly et al. [14] have experimentally studied adiabatic 
transitions of two electrons in a pair of tunnel- coupled 
QDs from a spin singlet S(2,0), with both electrons in 
the same QD and total spin 0, across an energy level 
anti-crossing to a spin triplet T+(l,l), with one elec- 
tron in each QD and total spin 1 (see Fig. [1]). Due to 
angular momentum conservation, the hyperfine-induced 
transition S-T+ involves an electron spin flip accompa- 
nied by a nuclear spin flop. In the experiment, the pro- 
cess was repeated many times with intermediate fast re- 
setting to the singlet state (200 ns). If the slow S-T+ 
transition was fully adiabatic and the nuclear spin po- 
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FIG. I: (color online) (a) Energy diagram for the relevant 
states of the double QD as a function of the bias e. (b) 
The hyperfine interaction allows for fiip-flop processes that 
open a gap Ahf at the S-T+ crossing. Driving the system 
slowly from S(2,0) can result either in an adiabatic (b) or a 
non-adiabatic (c) transition. Driving the system back very 
fast results in a non-adiabatic passage to two distinguishable 
charge states. 



larization sufficiently long-lived, this cycle should allow 
for complete nuclear polarization as the number of cycles 
becomes comparable to the number of nuclear spins in 
the QDs (typically about 10^ to 10^). However, a po- 
larization of only about 1% was achieved. Nevertheless, 
the coherence time measured in the S-Tq subspace was 



improved by a factor of up to 70 [1J| which can be at- 
tributed to the preparation of a suitable nuclear state. A 
previous theoretical model has been used to calculate the 
evolution of up to 36 nuclear spins per dot for an initial 
mixed state with fixed angular momentum per dot [3] ■ 

Here, we propose a theoretical model of nuclear state 
preparation taking into account all possible angular mo- 
menta, up to 200 nuclear spins per QD, and the possi- 
bility of partially adiabatic transitions. Such transitions 
can explain the smallness of the nuclear polarization gen- 
erated while allowing for a weak measurement ,16] of the 
nuclear state. The repeated weak measurements at each 
cycle lead to nuclear state narrowing and, as a conse- 
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quence, to a prolongation of the spin singlet-triplet co- 
herence. The main feature of our theoretical model is the 
description of the partially adiabatic S-T^ transitio n by 



the Landau-Zener-Stiickelberg (LZS) theory jlTI. llSl. Il9l|. 
Due to the high degeneracy of the nuclear spin state and 
the finite sweep tinie, we use suitable generalizations of 
the LZS model 

Model — We derive the Hamiltonian describing the par- 
tially adiabatic dynamics at the S-T+ anti-crossing. We 
start with the spin-preserving part Hq of the Hamilto- 
nian, with [Ho,Sz] = 0, describing the coupling between 
two electrons in a double QD in a magnetic field B, 



(1) 



where g* denotes the effective Lande g-factor (—0.44 for 
GaAs), fj,B the Bohr magneton, and the indices i ~ 1,2 
and a =t,i label the dot number and spin. The first 
term is the single-particle energy of the confined elec- 
trons, the second accounts for the Coulomb energy u of 
two electrons on the same QD, and the last for the elec- 
tron tunneling with strength r between the dots. 

The diagonalization of the first two terms of Eq. ^ 
leads to the relevant states of a double QD: the singlets 
S(0, 2), S(2, 0), S(l, 1) and triplets T±,o(l, 1), where (/, r) 
indicates the number of electrons in the {left, right) dot. 
The other states can be neglected as they have energies 
much higher than those considered here. The degeneracy 
of the singlets S(2,0) and S(l, 1) at e = ±w is lifted by 
the inter-dot tunneling, resulting in a splitting of \/2t. 
The energy levels as a function of the bias e = ei — £2 
are shown in Fig. [Ija). At the degeneracy e = e* of 
the singlet S(l,l) and the triplet T+(l,l), the hyper- 
finc interaction between the electron spins Si and the 
nuclear spins opens a splitting Ahf- The contact hy- 
perfine Hamiltonian is given by Hup = Si • h2 + S2 • h2 
where h; = X]fc=i ^i^i the Overhauser (effective nu- 
clear field) operator. The sum runs over the n nuclear 
spins in dot i, — Ufci^o |^'(rfc)P is the hyperfine cou- 
pling constant with the fc-th nucleus in dot i, with 'J'(rfe) 
the electron wave function, vq the volume of the unit 
cell and Vk the hyperfine coupling strength. Introducing 
Sf = Sf ± iSf and hf = ftf ± i/if , we write -Shf as 
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Since r 3> Ahf, we can use H — Hq + Hup to de- 
rive an effective Hamiltonian for the subspace spanned 
by {|S, ji,mi, j2,TO2)} and {|T+, ji, mi, ^2, m'a)}, where 
ji is the total nuclear angular momentum in dot i and 



nii its projection along B, 

H{e) ^Y.E,^{e)\qx){qx\ + \Y.{StK +Srht) (3) 
qx i 

where 9 = S,T, \x) = |ji, "t-i, ^2, "I2), and £^s,x(*^) = 
Esie), a.nd ET,^ie) = Et{s) + {hf + ^)x/^ + 9* (^bB. 

Method — With a time-dependent bias e — e{t), the 
Hamiltonian ([3]) is of the form H — Ho{t) + Hint, as 
the one studied by LZS to derive the staying and tran- 
sition probabilities Pa and Pna between two levels |1) 
and 1 2) driven through resonance between = — 00 and 
tf — -|-cx) by assuming that their energy difference is a 
linear function of time, A(i) = |i?s,x(0 ~ ^t.x(OI = 
at with the well-known result Pa = 1 ^ ^na = 1 ~ 
exp(— 27r|(l|i/HF|2)P/Q;S). Here, we study the effect that 
the transitions can have on the nuclear difference field 
8h^- = h\ — h\ and more precisely on its fluctuations 

cr^^) = \j {(5h^Y) — ((5/i^)^ that are responsible for the 

qubit dccoherence via T:^ = hjo^''^ [llii]. 

To study the evolution of cr^^-*, we compute the new 
state of the system after each cycle, which consists of a 
forward and a return sweep. The state of the system 
after the forward sweep results from the time evolution 
operator generated by ([3]); the resulting state is taken as 
the initial condition for the return sweep which is per- 
formed sufhciently fast to ensure a sudden parameter 
change without change in the state. As in the original 
LZS model, we assume a linear dependence A{t) — at 
throughout. Moreover, to make the problem treatable, 
we assume a constant hyperfine coupling Ak ~* A and 
the nuclear spins to be 1/2 (in reality the nuclear species 
of GaAs have spin 3/2). Before treating the case of many 
nuclear spins, we explain the main ideas by consider- 
ing the simple case of one nuclear spin in dot 1 and 
none in dot 2. The initial nuclear state is assumed to 
be an incoherent mixture of spin up and down; its den- 



sity matrix p^"^ has matrix elements — Sij/2 with 
hj =Tii- The initial mean value and standard devia- 



tion for Sh^ are {Sh^)o — and CTq = A/2. We now 
assume that after a cycle, a measurement of the electron 
spin (via charge) is performed to determine if a flip-flop 
has occurred. If a singlet is detected, the nuclear den- 



sity matrix becomes — e ^/2Ps, p^^^' — l/2Ps, 

and Pii"*'^ = P^V\ '^ ~ with 77 = nA^/ah and the 
probability to measure a singlet, Pg = (e^** + l)/2. In 
this new state, = ^(e^** - l)/2(e"'' + 1) and 

cr^""^ = (A/2)cosh"^(77/2). In the "fast" limit a ^ 00, 
we find pn"^'^ — pn"' therefore the variance is unchanged, 
a{ = ' . In the "slow" limit a ^ 0, we have 
P^ll'^ ~ ^ while all other elements of pn'^'^ vanish, such 

that crj^' — 0, describing a strong (projective) measure- 
ment of the nuclear spin. For < a < 00, the detection 
of a singlet induces a weak measurement which decreases 
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the fluctuations, uj^'' < (Tq^^ On the other hand, if a 
triplet is measured, we have p\^^'"^ = P^i'"^ ~ p'l^'"^ — 

(1) T 

and pj^j^ =1 independently of a; in this case the nu- 
clear spin is projected on the down state. We conclude 
that it is possible to reduce the fluctuations of the nuclear 
spins without fully polarizing them. Below, we show that 
the same mechanism works also for a system with many 
nuclear spins. 

For n ^ 1, the nuclear states are highly degener- 
ate and the LZS propagator from the simple case can- 
not be used anymore. An elegant solution to derive the 
LZS propagator for degenerate systems consists in ap- 
plying the unitary Morris-Shore (MS) transformation to 
the LZS Hamiltonian to reduce the dynamics into sets of 
decoupled single states and independent two-level sys- 
tems [20]. Because i?HF exclusively couples states of 
the form \S, ji,mi, j2,m2) = |0) to the degenerate states 
|T+,ji,mi-l,j2,nT-2> = |1) and |T+, ji,TOi, J2,TO2-1) = 
1 2), the LZS Hamiltonian can be brought into a block di- 
agonal form. The MS basis is found by diagonalizing 
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where Vij = {i\Hint\j), at = y/jt{ji + 1) - rn^{m^ - 1). 
The eigenstates are |0) and |1') = (ai|l) -|-a2|2)) /gq as- 
sociated with the eigenvalue = A^/Suq, and |2') 
(ai|l) — a2|2)) /ao with eigenvalue 0, where Oq = 
The states with the same eigenvalues are coupled with 
strength A. The state |2') is a "dark state", as it does 
not couple to other states. In the subspace spanned by |0) 
and the time-dependent Schrodinger equation with 
the initial state |0) at time <i can be solved and thus the 
LZS propagator elements can be calculated :21j]. In 
order to express the solution in the original basis, we 
perform the inverse MS transformation to find the ma- 
trix elements Uoi = {ai/ao)U^9. We also account for 
finite time propagation and thus avoid the unphysical 
situations of infinite energy that arises for couplings that 
do not vanish when t — > ±oo, and infinite detuning as we 
assume A{t) = at. The finite-time solution also allows 
us to model the situation in the experiment where 
mixing between To(l, 1) and the S(l, 1) states must be 
avoided. 

At typical operating temperatures and fields, where 
ksT ^ ^nMn-B, the initial nuclear density matrix 
can be assumed to be diagonal, p^ '' — 
with a uniform distribution of states x- 120 mK 

and 100 mT, we have k^T ~ 10~^eV and g^pj^B ~ 
10~^ eV. The joint probability p{x) can be factorized 
into TOi)p'(j2, TO2) since the dots are initially in- 

dependent, with p'{ji,mi) = g{ji)f{mi\ji), where g{ji) 
is the probability for total nuclear angular momentum 
ji and f{m,\ji) = [e{j^ + m,) - 0{j, - m,)]/{2j, + 1) 
is the equally distributed conditional probability of hav- 



ing a magnetization uii given ji, and is the Heaviside 
function with 9{0) — 1. The probability distribution 
g{ji) — G{ji)/J2ji ^(■?») found by counting how many 
times G{ji) an irreducible representation of dimension 
2ji + 1 occurs. 

After the forward sweep of the fc-th cycle, the state is 
^(fc+i) ^ Up^^W\ with p(o) = |S)(S| p'i\ The back 
sweep will act as a measurement of the final configura- 
tion of the electronic system. After the LZS transition 
the charge configuration of the system is (1,1) indepen- 
dent of the spin, but after the back sweep it is a superpo- 
sition of T+(l, 1) and S(2,0) which has a relatively fast 
decay time t '--^ 1 ns [l^l , such that the electronic system 
will evolve with probability Ps to the singlet S(2, 0) and 
with probability Pt to the triplet T+(l,l) after a time 
T. This provides a way to determine if the system has 
evolved adiabatically or not during the forward sweep. 



We write p^'^+i) = PsPs""^^ 



PtPx^^' to describe the 



mixture of the ensembles that have evolved adiabatically 
and non-adiabatically in the forward sweep. After a time 
T, this state will collapse either to or with 

probability Ps and Pp, respectively, with [q ~ S, T), 



P„ = Tr 



MqUp^'^'lUHll 



(5) 



where Mq = \q){q\ is the projection operator describing 
a strong measurement in the charge sector. If a singlet 
{q = S) or triplet {q = T) is detected after the (fc -I- l)-th 
sweep, we update p^ according to 



(6) 



In the case of a triplet {q = T), we must also take into 
account the nuclear spin flop. 



(fc+i) 

p 

•31323132 
p^'+'^ 

f^Jl32rnim2 



^0, 



(fc+l) (fe-l-1) 

/^Jij2,mi + l,m2 Pjij2'mi,m2 + l 



^3l32mim2 yPjij2-mi + l,m2 ~^ Pjij2mi,m2 + 1 



(7) 

After a flip-flop process, the crossing point e* moves 
slightly because of the change in the total magnetic fleld 
B -I- (hi -I- h2)/g*pB and affects the initial and flnal times 
of the LZS transitions. With the convention that i = at 
the initial anti-crossing point, the time shift is U ti+St, 
t{ + St, with St = A/a. 

We have performed a Monte Carlo simulation where 
ti — —30 ns and tf = 20 ns, such that the duration of the 
flrst sweep is the same as in [13] • Our choice for a satisfies 
the conditions A <C ati and A <C atf <^ g* p^B that hold 
for a system far from resonance at the beginning and at 
the end of the first cycles, and with a final difference in 
energy smaller than the Zeeman splitting. The order of 
magnitude of A is determined by looking at the case of the 
most likely ji, which is ^ ^Jn, with = to consider 
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FIG. 2: (color online) (a) Averaged Overhauser field fluctua- 
tions cr'^-' and electron spin decoherence T2 = h/cr^^'^ (upper 
inset) as a function of the number of performed cycles, (b) 
Singlet return probability Ps as a function of the number of 
performed cycles. Error bars are smaller than line thickness, 
(c) Initial (black) and final (orange) averaged probability dis- 
tribution of the nuclear spin eigenstates of Sh'' . Repeated 
cycles narrow the distribution, showing a reduction of o"*-^'. 



cycling combined with a spin-to-charge conversion to the 
(1,1), respectively (2,0), configuration induces a weak 
measurement on the nuclear state which strongly con- 
tributes in the suppression of the nuclear fluctuations. 
An experimental confirmation of the predicted weak mea- 
surement would be provided by measuring Pg « 1 and 
at the same time an increase in Tj* as a function of the 
number of cycles. 

Using a Monte-Carlo algorithm for n ~ 200 nuclear 
spins per QD, we find an enhancement of the electron 
spin coherence time Tj* = h/a^^^ by a factor of ~ 1.34 
after 300 cycles. While the initial al oc y/n, the fi- 

(z) 

nal cTj is mainly determined by the preparation mech- 
anism and is approximatively independent of n. There- 
fore, we expect T2i/T2[ = jaY^ oc \pn^ and we es- 
timate T* JTI^ w 94 for n w 10^. This scaling stops 
when T2* w T2 or before if dipolar interaction effects are 
taken into account. Future calculations for more spins 
will allow more direct comparison with experiments. Our 
results also apply for QDs with fewer nuclear spins e.g. 
in Si, ZnO, carbon nanotubes, or graphene. 
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the strongest coupling in that subspace. The number 
of spins is currently limited to n = 200 per dot by the 
computational power at our disposal. This limitation 
has a side effect on the choice of the hyperfine coupling 
constant. For large spin systems, one has A = ^ with 
A = 90 /ieV for GaAs, but for dilute systems A can be 
crucially smaller, therefore we choose A = 9neV 24 1. 
Thus, A ~ 10~® eV so that atf must be of the order of 
10"'^ eV since |5*^bS| ~ 2.3 ^eV. This implies that a 
must be between 10 and 50eVs~^, which also satisfies 
the condition on ati. 

Results — In Fig. (Ha), we plot the evolution of a^''\ 
averaged over 160 runs, as a function of the number of 
cycles for a = lleVs"^. Since the reduction of cr*^^^ 
persists in the average, a read-out of the charge state 
after each cycle is not required. When the polarization 
reaches about 8% (generally at « 1/ ^Jnj2) the nuclear 
spins approach a "dark state" preventing further polar- 
ization. We plot in Fig. [2l^b) the average probability Pg 
of measuring a singlet as a function of the number of cy- 
cles. Ps can be measured in the S - T+, as it was done 
for S - To 0. 

Conclusion — We have developed a model to explain 
the increase of T2 in double QDs via the tuning of two 
electrons across an energy level anti-crossing between 
S(2, 0) and T4.(l, 1). Our model is based on the possibil- 
ity of partially adiabatic transitions which are described 
by a generalized LZS theory. We have shown that the 
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